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MEASUREMENT  OF  ENERCY  FLUX  DENSITY  DISTRIBUTION 


IN  THE  FOCUS  OF  AN  ARC  IMAGE  FURNACE 


Rodney  B.  Beyer,  Leonard  McCulley,  and  Marjorie  W.  Evans 


SUMMARY 

A  calorimeter  for  mapplni;  flux  density  In  a  focal  volume 
la  described.  The  deslgrn  Is  based  on  an  earlier  deslfrn  by 
Broldo  and  Wllloufthby.  It  differs  chiefly  In  two  ways.  First 
It  Is  much  smaller  In  order  to  permit  Its  placement  Inside  a 
test  specimen  chamber.  Second  the  deslftn  provides  for  an 
aperture-containing  shield  which  Is  placed  In  front  of  the 
receiving  surface  calorimeter  In  order  to  limit  the  area  of 
beam  which  Is  Intercepted.  The  decay  constant  Is  2  x  10'^ 
sec~‘.  Flux  densities  ranging  from  5  to  120  cal/cm*-sec 
can  be  conveniently  measured  over  Intercepted  beam  areas  s 
of  from  0.03  to  0.2  cm*. 

With  this  calorimeter  the  flux  density  distribution  In 
the  focal  volume  of  a  double  ellipsoidal  arc  Image  furnace 
system  was  determined. 
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Contract  Nonr-341S(00) . 


1 


The  arc  image  furnace  Is  a  useful  device  for  supplying  energy  at 
controllable  flux  density  levels  and  for  controllable  times. For 
example,  in  work  soon  to  be  published  it  has  been  used  in  studying  the 
initiation  of  dellagratlon  waves  at  solid  surfaces.  The  radiation  beam 
from  the  douhle-el  11 psol dni  carbon  arc  image  furnace  used  in  this  work, 
shown  schematically  in  Fig.  1,  converges  toward  a  1-cm  diameter  focus  in 
a  cone  with  vertex  angle  of  120  degrees.  Ad.justable  shutters  (not  shown 
in  Fig,  l)  placed  near  the  sample  reflector  allow  the  magnitude  of  flux 
density  at  the  focus  to  be  ad.Justed.  The  focus  is  an  ellipsoidal  volume  in 
which  the  flux  density  varies.  For  the  furnace  to  be  a  useful  instrument 
in  a  quantitative  study,  one  must  know  the  flux  density  d-*  stribution 
throughout  this  focal  volume  and  especially  the  maximum  flux  density. 

The  purpose  of  this  work  was  to  design  a  suitable  calibrating  calorimeter, 
and  with  it  to  map  the  flux  density  distribution. 

1 .  Design 

The  calorimeter  design  was  an  adaptation  of  one  by  Broldo  and 
Willoughby.®  The  receiver  was  constructed  from  a  copper  disc,  0.10  cm 
thick  and  approximately  0.4  cm  in  radius.  This  was  pressed  into  a 
spherical  segment  with  an  Inner  radius  of  0.55  cm  and  an  inner  height  of 
0.15  cm.  The  Inner  surface  was  blackened  with  electro-deposited  platinum 
black.  A  plane  silver  foil  shield,  0.0127  cm  thick,  with  an  aperture 
which  was  varied  between  0.025  and  0.508  cm  in  diameter,  was  placed  in 
front  of  the  concave  surface  to  limit  the  area  of  the  beam  intercepted 
by  the  receiver.  The  shield  was  spaced  0.025  cm  from  the  rim  of  the 
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receiver,  with  the  line  connecting  the  center  of  the  aperture  and  the 
center  of  the  receiver  parallel  to  the  optical  axis  of  the  system.  A 
5>mll  copper-constantan  thermocouple  was  silver-soldered  to  the  back  side 
of  the  receiver  near  one  edge  to  measure  the  temperature  rise.  The 
receiver  was  mounted  on  Bakellte  posts  for  thermal  Insulation  from  the 
holder  and  the  plate  containing  the  aperture.  A  sketch  of  the  receiver 
and  the  shield  Is  shown  in  Fig.  2.  The  physical  constants  of  the 
receiver  were 

a  B  absorptivity  «  0.97 

p  <=  density  b  8.92  gm/cm’ 

c  =  specific  heat  =  0.093  cal/®C-gra 

X  =  thermal  conductivity  =  0.91  cal/cm-®C-8ec 

1C  =  X/p  c  ■  thermal  diffusivity  »  1.095  cm*/sec 

£  thickness  0.10  cm 

m  m  mass  =  0.406  gm 


2.  Theory 


If  heat  losses  can  be  neglected,  if  the  gradient  across  the  thick¬ 
ness  is  small,  and  if  the  difference  between  the  temperature  measured 
at  the  edge  and  the  average  temperature  is  small,  the  average  flux  density 

I  at  the  aperture  over  an  aperture  area  A  is  given  over  any  time 
a  a 

interval  At  by 


I 

a 


m  c 
a  A 


A6 


AV 


A* 


(1) 


a  \  I 

where  is  the  increase  of  average  temperature  of  the  receiver,  and 

the  temperature  may  be  measured  anywhere  on  the  calorimeter.  The  flux 
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density  versus  time  curve  of  a  pulse  of  energy  is  obtained  by  differen¬ 
tiating  the  temperature-time  curve  of  the  calorimeter. 


Following  Broido  and  Willoughby  the  characteristics  of  the  calorim¬ 
eter  are  defined  in  terms  of  (a)  the  decay  constant  A  (sec“‘), 

h 


A  pc6 
c 


(2) 


where  h  is  the  heat  loss  coefficient  of  the  receiver  in  cal/sec-*’c  and 

A  is  the  calo: 
c 

stant  T  (  sec) , 


A^  is  the  calorimeter  surface  area^  and  (b)  a  differentiated  time  con- 


1.7fi* 

Kn* 


(3) 


They  showed  that  if  At/2  <  <  1  ,  heat  losses  can  be  neglected.  Again 
following  Broido  and  Willoughby  h/A^  was  estimated  by  assuming  con¬ 
ductive  heat  losses  to  be  negligible  and  convective  heat  losses  to  be 
equal  to  radiative  losses,  as  calculated  from  the  Stefan-Boltzmann  law. 

In  our  measurements  the  calorimeter  temperature  rise  was  varied  between 
10  and  50®C  and  was  normally  20®C.  For  a  rise  of  20®  from  room  tempera¬ 
ture  of  27®C,  A  =  2  X  10"*  sec"*.  Thus  for  exposures  to  pulses  even  as 
long  as  our  maximum  of  10  seconds,  the  heat  losses  may  be  neglected. 
Oscillographic  traces  of  the  calorimeter  temperature  showed  the  cooling 
curve  slopes  to  be  negligible,  in  agreement  with  this  calculation.  Broido 
and  Willoughby  further  showed  that  for  t  small  with  respect  to  the  expo¬ 
sure  time,  the  temperature  gradient  across  the  thickness  6  is  small, 
and  the  measure  of  the  rate  of  change  of  back  surface  temperature  for  a 
uniformly  exposed  receiver  Is  a  valid  measure  of  the  rate  of  change  of 
the  average  temperature.  For  this  receiver  t  Is  0,0016  sec,  while  flux 
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measurements  were  made  with  exposure  times  from  0.030  to  10  sec.  Thus 
^/At  can  be  measured  at  the  back  surface.  AO/At  was  constant  over  the 
entire  range  of  temperature  rise  and  exposure  time  for  each  of  the  aper- 
tures;  showing  that  heat  conduction  or  convection  from  neighboring  parts 
of  the  apparatus  was  negligible. 

The  received  flux  density  in  our  case  was  not  distributed  uniformly 
over  the  receiver  surface.  This  was  because  the  flux  density  varied  within 
the  focal  volume  and,  more  importantly,  the  area  receiving  the  flux  idiich 
passed  through  the  aperture  was  in  general  smaller  than  the  total  surface 
area  of  the  receiver.  Thus  it  was  necessary  to  show  that  the  temperature 
measured  by  the  thermocouple  at  the  edge  of  the  calorimeter,  0^ ,  was 
to  good  approximation  the  average  temperature  0^^  .  This  can  be  tested 
by  considering  the  problem  of  a  disk  of  diameter  R|  and  thickness  6  ) 
heated  over  a  circular  area  of  radius  Rj  <  R)  . 

The  flow  of  heat  in  such  a  circular  disk  is  represented  by  the  fol¬ 
lowing  differential  equation  where  it  is  assumed  that  the  disk  is  suffi¬ 
ciently  thin  that  only  a  radial  flow  exists; 

®  ^R  1 

®t  =  77#  ^  ^ 

with  ■  0  ■  0  at  t  «  0 

9^  =  0  at  r  *  0,  Rj  , 

{1  X  >  0 

0  X  <  0 

Here  t  and  r  are  time  and  radial  independent  variables  and  the  sub¬ 
scripts  t  and  r  indicate  partial  derivatives,  t  is  exposure  time 
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is  the  flux  density  at  the  calorimeter  over  the  circle  of  radius 
R,  •  E^juatlon  (4)  is  solved  by  reducing  it  to  two  partial  differential 
equations : 


9^  ■  +  Z  9j 


rr  r  r 


for  Rj  <  r  <  R, 


(5) 


with 

e 

«  0 

at  t  •  0 

9 

r 

=  0 

at  r  •  R, 

and 

a 

P 

H(T-t)  +  ^(0^^  +  -  for  0<r<R, 

(6) 

with 

0 

>  0 

at  t  •  0 

0 

r 

•  0 

at  r  •  0  , 

On  taking  the  Laplace  transform 

e(a, 

r) 

ao 

«  I  e  0(  t,  r)  dt 

a)  0 

(7) 

of  Equations  (s)  and  (6)  along  with  their  boundary  conditions,  one 

obtains 

the 

following  ordinary  differential  equations  in  the  variable  r  : 

s9 

e 

7  '  R2  <  T  <  R, 

rr  r  r 

(8) 

with 

0 

r 

■  0  at  r  ■  Rj  , 

and 

a 

S0  *=  — 

P 

c  6 

^  ^  +  ic(9  +  -  0  )  ,  0  <  r  <  Rg 

s  rr  r  r  '  -  ' 

(9) 

with 

•  0  at  r  ■  0  , 

«  I_  ,  _  -ST 

Define  q*  ■  -  and  Of  m  - *— r  - - - 

tf  1C  p  c  6  s 

Equations  (8)  and  (9)  are 


•  'Rien  the  solutions*  of 
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R,  <  r  <  R,  (10) 


^  +  C,  lo  C  qr)  ,  0  <  r  <  Rj 


(11) 


where  C.  and  C.  are  arbitrary  constants  and  I  and  K  are  the 
*  *  n  n 

solutions  of  the  modified  Bessel  equation.  The  constants,  Cj  and  Cj, 
are  determined  by  applying  the  further  conditions  that  0  ,  in  ( 10) 

and  ( ll)  must  agree  identically  at  r  e  r,  ,  Equations  ( 10)  and  (ll)  become 

q  ll  i  q  Kt; 


R,  <  r  <  R, 


(12) 


and 


^  Y^lqR-J 

0  <  r  <  R*  .  (13) 

Set  r  ■  Ri  in  Equation  ( 12)  and  r  ■  0  in  Equation  ( 13) .  Then 

(12a) 


^  ^  (q 

q»  R,  I,  (q  R.) 


at  the  edge 


and 


0  ,  ®  +  a_5i  ["i  (q  R,)  /  ‘I  -  Kj  (q  Rj)!  at  the  center.  ( 13a) 

q  <>  L  ll  Ri^  J 

Laplace  inversion*  of  Equations  ( 12a)  and  ( 13a)  gives  Integrals  which 
may  be  evaluated  using  Cauchy's  theorem  of  residues.*  Let  R  •  R2/Ri  • 
Define 


Cp  (t,  R,)  « 


li 

p  c  6 


R*  +  R»t  -  2  RSj. 


and 


-  7^  [- 1?^ 


(14) 


(15) 
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where 


and 


"  i£l  K  o-i*  ^  Ri  Jo  (ffj  Ri) 


;  -Kai^t  J,  (o'  Rj) 

y t)  .  yt)  +  "ToTiTRj 


tti  R,,»‘-/r+l 


J 


H,f  kD^-lzir),  7  UV  1 

r-O  rj  ( r  +  1)1  \ifel  wp"! 


;  ^ 


and  Jo,  J,  are  the  ordinary  Bessel  functions.  The  are  the  non¬ 


negative  solutions  of  J,(o  R,)  •=  0.  Then 

^  (t,  Rj)  ■  CD  (t,  Rj)  -  H  (t-x)  CP  (t-T,  Rj) 

(t,  0)  «  CD  (t,  0)  -  H  ( t-r)  CD  (t-T,  O) 

Assume  t  -  t  <  0,  Then  Equations  ( 16)  and  (l7)  become 


0  ( t,  R,)  ■  CD  (t,  R,)  s  at  the  edge 

A  ( t,  0)  «=  CD  ( t,  0)  at  the  center  . 

At  time  t,  the  average  temperature  of  the  disk  is  given  by 


a  I  R»  t 
p  c  S 


(16) 

(17) 

(16a) 

(I7a) 

(18) 


The  ratio  of  the  temperature  at  the  edge  of  the  disk  to  the  average  tem¬ 
perature  at  any  time,  t,  is  found  by  combining  Equations  ( 14) ,  (l6a), 
and  ( 18)  to  give 

A, 
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1_ 

AV 


1 

8  V  t 


2Ej.(t) 

Rt 


(19) 


Note  that  lim  l!_(t)  ••  0,  It  is  therefore  possible  to  find  a  t,, 
t-=«»  E 

such  that  9jj  y  as  closely  as  desired  for  all  t  >  to  where 
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9. 


p  >  _  p  s 
1  _  “j - 5a- 


(20) 


R  ■  8  If  t 

The  arc  length  along  a  geodesic  line  across  the  calorimeter  surface  is 

0.828  cm.  For  an  aperture  diameter  of  0.25  cm  the  arc  length  subtended 

by  the  beam  is  0.604  cm.  In  approximation  take  Rj  «  0.414  cm  and 

R,  “  0.302  cm.  For  a  typical  case  with  this  aperture  =  20  cal/cm’-sec, 

and  At  >  1.17  seconds  for  ^  e  20°C.  From  Equation  (20)  A  s  0.992 

R 

for  If  c  1.095  cm^/sec.  Considerations  of  the  asymptotic  forms  for  the 

Bessel  functions  and  Inspection  of  the  leading  terms  allow  one  to  bound 

the  series  ^_( t)  so  that 
E 


9-  -  9„ 

R  R 


BEJt) 

— = -  <  e”‘®® 

t  R 


at  t  -  1.17  sec. 


Thus  the  error  in  approximating  6  by  6  in  this  case  is  negligible 
and  the  measurement  of  temperature  at  the  edge  of  the  calorimeter  gives 
a  good  approximation  to  the  average  temperature. 

In  summary  the  analysis  shows  that  the  average  flux  density  at  the 
plane  of  the  aperture  and  over  the  aperture  area  may  be  calculated  from 
Equation  (l),  where  is  to  good  approximation  the  rate  of  tem¬ 
perature  rise  at  the  thermocouple. 


3.  Flux  Density  Distribution  in  the  Focal  Volume 

Flux  density  in  the  focal  volume  was  mapped  by  traversing  along 
axial,  horizontal,  and  vertical  axes  with  constant  apertures  in  the  range 
0.025  to  0.508  cm  in  diameter.  In  further  experiments  the  calorimeter  was 
centered  at  a  selected  point,  and  the  flux  density  was  measured  for  dif¬ 
ferent  apertures  over  the  same  range.  Figure  3  shows  flux  density 
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normalized  by  maximum  flux  density,  as  a  function  of  vertical,  horizontal, 
and  axial  distance  from  the  point  of  maximum  flux.  The  aperture  was 
0.12'^  cm.  For  the  horizontal  traverse,  points  for  an  aperture  of  0.064  cm 
are  also  shown;  the  shape  of  the  surfaces  of  constant  flux  density  was  the 
same  for  each  aperture.  The  dispersion  is  believed  due  chiefly  to  wandering 
of  the  arc  and  hence  of  the  image  over  the  aperture. 

A  typical  curve  of  measured  flux  density  versus  aperture  diameter, 
with  the  calorimeter  in  this  case  centered  at  the  point  of  maximum  flux 
density,  is  given  in  Fig.  4.  For  aperture  diameters  smaller  than  0.250  cm 
the  apparent  flux  density  rises  sharply.  Since  for  all  apertures  the  flux 
density  drops  off  as  one  moves  away  from  the  maximum  at  the  focal  volume 
center,  the  increase  in  apparent  flux  density  with  decrease  in  aperture 
at  given  calorimeter  position  must  be  related  to  aperture  thickness.  At 
an  aperture  of  0.250  cm  the  aperture  thickness  begins  to  be  a  significant 
fraction  of  the  radius,  the  ratio  being  0.1.  A  thickness-related  increase 
in  apparent  flux  density  can  arise  from  reflections  from  the  walls  of  an 
aperture  placed  in  a  focal  volume.  From  geometric  optics  it  is  clear  that 
proper  placing  of  a  thick  aperture  within  the  focal  volume  can  substantially 
increase  the  flux  passing  through  the  aperture  above  the  amount  which  would 
pass  if  the  shield  were  infinitely  thin.  Auxiliary  experiments  confirmed 
that  apparent  flux  density  was  dependent  on  foil  thickness  at  small 
apertures. 

The  data,  which  include  Figs.  3  and  4  as  well  as  traverses  with  larger 
apertures,  are  consistent  under  the  assumption  that  effects  due  to  reflec¬ 
tion  from  the  wall  of  the  aperture  are  negligible  at  apertures  of  0.25  cm 
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normalized  by  maximum  flux  density,  as  a  function  of  vertical,  horizontal, 
and  axial  distance  from  the  point  of  maximum  flux.  The  aperture  was 
0.12'7  cm.  For  the  horizontal  traverse,  points  for  an  aperture  of  0.064  cm 
are  also  shown;  the  shape  of  the  surfaces  of  constant  flux  density  was  the 
same  for  each  aperture.  The  dispersion  is  believed  due  chiefly  to  wandering 
of  the  arc  and  hence  of  the  image  over  the  aperture. 

A  typical  curve  of  measured  flux  density  versus  aperture  diameter, 
with  the  calorimeter  In  this  case  centered  at  the  point  of  maximum  flux 
density.  Is  given  in  Fig.  4.  For  aperture  diameters  smaller  than  0.250  cm 
the  apparent  flux  density  rises  sharply.  Since  for  all  apertures  the  flux 
density  drops  off  as  one  moves  away  from  the  maximum  at  the  focal  volume 
center,  the  Increase  In  apparent  flux  density  with  decrease  in  aperture 
at  given  calorimeter  position  must  be  related  to  aperture  thickness.  At 
an  aperture  of  0.250  cm  the  aperture  thickness  begins  to  be  a  significant 
fraction  of  the  radius,  the  ratio  being  0.1.  A  thickness-related  Increase 
In  apparent  flux  density  can  arise  from  reflections  from  the  walls  of  an 
aperture  placed  In  a  focal  volume.  From  geometric  optics  It  Is  clear  that 
proper  placing  of  a  thick  aperture  within  the  focal  volume  can  substantially 
Increase  the  flux  passing  through  the  aperture  above  the  amount  which  would 
pass  If  the  shield  were  infinitely  thin.  Auxiliary  experiments  confirmed 
that  apparent  flux  density  was  dependent  on  foil  thickness  at  small 
apertures. 

The  data,  which  Include  Figs.  3  and  4  as  well  as  traverses  with  larger 
apertures,  are  consistent  under  the  assumption  that  effects  due  to  reflec¬ 
tion  from  the  wall  of  the  aperture  are  negligible  at  apertures  of  0.25  cm 
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and  greater .  They  depict  an  image  with  roughly  ellipsoidal  surfaces  of 
constant  flux,  the  maximum  flux  density  being  at  the  center .  The  flux 
density  falls  off  to  97-99^  of  maximum  at  0.05  cm  from  the  center,  and 
to  88-95^  of  maximum  at  0.125  cm  from  center.  The  maximum  value  is  5^ 
higher  than  that  measured  by  the  0.250-cm  aperture  centered  at  the  focal 
center . 
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FIGURE  TITLES 


1.  Sketch  of  double  ellipsoidal  arc  Image  furnace  system. 

2.  Calorimeter  for  calibrating  radiant  flux  density. 

3.  Flux  density  divided  by  maximum  flux  density  versus 
distance  from  center  of  focal  volume. 

4.  Apparent  flux  density  at  center  of  focal  volume  versus  aperture 
diameter. 
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FIG.  1  SKETCH  OF  DOUBLE  ELLIPSOIDAL  ARC  IMAGE  FURNACE  SYSTEM 


FIG.  2  CALORIMETER  FOR  CALIBRATING  RADIANT  FLUX  DENSITY 
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FIG.  3  FLUX  DENSITY  DIVIDED  BY  MAXIMUM  FLUX  DENSITY  VERSUS 
DISTANCE  FROM  CENTER  OF  FOCAL  VOLUME 


aperture  diameter  —  cm 
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FIG.  4  APPARENT  FLUX  DENSITY  AT  CENTER  OF  FOCAL  VOLUME  VERSUS 
APERTURE  DIAMETER 


